Two— step condensation of the charged Bose gas 
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The condensation of the spinless ideal charged Bose gas in the presence of a magnetic field is 
revisited. The conventional approach is extended to include the macroscopic occupation of excited 
kinetic states lying in the lowest Landau level, which plays an essential role in the case of large 
magnetic fields. In that limit, signatures of two diffuse phase transitions (crossovers) appear in the 
specific heat. In particular, at temperatures lower than the cyclotron frequency, the system behaves 
as an effectively one-dimensional free boson system, with the specific heat equal to ^Nks and a 
gradual condensation at lower temperatures. 



I. INTRODUCTION 

The equilibrium behavior of a charged Bose gas (CBG) 
in the presence of a magnetic field is a fundamental prob- 
lem in quantum statistical mechanics. It was first ad- 
dressed by Schafroth [11] at a time where the BCS the- 
ory of superconductivity had not yet been discovered. 
Although it was soon recognized that a model of non- 
interacting bosons had little predictive power for super- 
conductivity, the CBG problem has continued to attract 
considerable attention because of its appealing simplicity 
and as a model to simulate astrophysical scenarios where 
charged particles are subject to extremely large magnetic 
fields, or neutral particles move in a fast rotating back- 
ground. It is also interesting as a preliminary step in 
the study of the condensation of interacting bosons in a 
magnetic field or under rotation. 

For the non-interacting case, Schafroth noticed that, 
in the presence of a magnetic field, Bose-Einstein con- 
densation (BEG) no longer occurs in a strict sense, since 
at low enough temperatures, the system behaves as one- 
dimensional (ID) and a non-interacting ID boson system 
is known to exhibit not a sharp but a gradual (diffuse) 
BEG transition. Later, May [J Q investigated the con- 
densation in the non-relativistic case and showed that a 
sharp transition occurs for dimension > 5. Daicic and 
coworkers [3, Hj extended May's findings to the relativis- 
tic high-temperature case. Using a different definition of 
BEG, Toms proved that BEG in the presence of a 
uniform magnetic field docs not occur in any number of 
spatial dimensions and Elmfors and coworkers consid- 
ered that, in the three-dimensional (3D) case, although a 
true condensate is not formed, the Landau ground state 
can accommodate a large boson number. In this spirit, 
Rojas [m found that BEG may occur in the pres- 
ence of a homogeneous magnetic field, but there is no 
critical temperature at which condensation starts, the 
phase transition being diffuse. Diffuse phase transitions 
are those not having a definite critical temperature, but 
an interval of temperatures along which the transition 
occurs gradually. The concept was already introduced 
in the study of phase transitions which occur in certain 
ferroelectric materials In this paper we treat the 

terms diffuse transition and crossover as synonyms. The 



notion of BEG of a GBG in the presence of a magnetic 
field in 3D continued to be studied by Toms and cowork- 
ers [ll,[l3], who showed that, although there is no BEG 
in the sense of a sharp phase transition, the specific heat 
exhibits a clear maximum that can be used to define a 
critical temperature. Using that definition, the authors 
inferred that the critical temperature increases with the 
magnetic field reaching the usual value for the 3D 
free Bose gas when no magnetic field is present. The ex- 
tension of these studies to the case of vector bosons in 
strong magnetic fields and at high temperatures was per- 
formed by Khalilov and coworkers p5l - [l7l |. Specifically, 
the condensation and effective magnetization of a charged 
vector boson gas were studied, showing that there is no 
true BEG as well, although a significant amount of bosons 
can accumulate in the ground state at low temperatures. 
More recently, the magnetic properties of charged spin-l 
Bose gases in an external magnetic field have been re- 
visited, with focus on the c omp etition between diamag- 
netism and paramagnetism [l8|. 

In the particular case of a large magnetic field, the 
main feature of a GBG is the coexistence of highly de- 
generate Landau levels and a fine structure of levels with 
nonzero momentum parallel to the magnetic field. We 
propose that a GBG under a large magnetic field can also 
be used as a crude model to understand basic features 
of systems where a fine sublevel structure coexists with 
sectors characterized by highly discrete quantum num- 
bers. Specifically, we have in mind the condensation of 
bosonic molecules, whose level structure is hierarchical, 
with their rotational states structured into vibrational 
sectors which in turn can be grouped into largely spaced 
electronic levels. 

Thus we are motivated by the study of molecular 
BEGs, where a complex molecular level structure may 
give rise to multi-step BEG. Molecular condensates are 
expected to display a wealth of fundamental phenom- 
ena. For instance, it has been proposed that molecular 
condensates may permit the experimental study of low- 
energy parity violation [l9|. Specifically, the fact that 
recently both homonuclear |20| . heteronuclear [2l|, [22j 
ultracold molecules and heteronuclear molecular ions 
[23I [2^ have been produced in the rovibrational ground 
state, suggests that molecular BEG may not lie too far in 
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the future. The conclusion is that the study of muhi-step 
condensation of the CBG is interesting not only from an 
academic point of view, but it may provide a qualitative 
understanding of the condensation of molecules with a 
hierarchical level structure. 

We consider a wide range of magnetic fields because we 
are interested in understanding the general trends. This 
means that in our analysis we include magnetic fields so 
large that would be unrealistic in some contexts [25l - l27| . 
Comparing our results with the previous related litera- 
ture on the CBG (see, for example, Refs. [H and [13), 
we notice the existence of a term that has so far been 
neglected. This is the occupation of states in the low- 
est Landau level but with arbitrary kinetic energy in the 
direction parallel to the magnetic field. This term is par- 
ticularly important at large magnetic fields and in fact 
is responsible for the two-step condensation we refer to 
in the title and which we shall discuss in detail. We 
advance here that, as the temperature descends, a first 
step is defined by the onset of a preferential occupation 
of states in the lowest Landau level but without any of 
them absorbing a large fraction of the bosons. Then a 
plateau is reached where the system behaves as effectively 
one-dimensional. The second condensation step occurs 
at even lower temperatures, when the true one-particle 
ground state becomes host to a macroscopic fraction of 
the total boson number. By the time this occurs, the sys- 
tems behaves as one-dimensional, which translates into 
a diffuse Bose-Einstein condensation, as opposed to the 
conventional, sharp phase transition which is character- 
istic of three dimensions. Interestingly, these two dis- 
tinct, both gradual condensation steps have been noted 
by Ketterle and coworkers [2^,[2^ in a different but some- 
what analogous context, namely, that of condensation in 
a strongly anisotropic harmonic trap. As compared to 
a CBG in a large magnetic field, the role of the Landau 
levels is played there by the transverse subbands char- 
acterizing the motion in the most confined direction. In 
Refs. and [2^ the above mentioned term of the oc- 
cupation of the states with the lowest value of the most 
discrete quantum number, was correctly included, which 
resulted in the prediction of a two-step condensation. 



rection, and both n and run over natural numbers 
(we assume hard wall boundary conditions at z = and 
z ^ L). The index n is said to characterize the Landau 
level. 

The area A is assumed to be large enough for the role 
of edge states in the x — y plane to be negligible. Specif- 
ically, P <C where P ~ hc/qB is the magnetic length 
squared. 

We start by writing the total number of particles as 
{h=l) 
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where /3 = (ksT)^^ , d = A/P is the Landau degeneracy, 
and z = cxp{/3fi) is the fugacity, with fi the chemical po- 
tential [which absorbs the 1/2 term in ([1])]. We partition 
the total particle number into four different groups: 



N = No+Nk+N„,+ Nrak ■ 



(3) 



Nq is the number of particles in the one-particle ground 

state {n ~ Uz — 0), 



Nn 



(4) 



Nk is the number of particles in the ground state of the 
magnetic level but in excited kinetic levels (n = and 
riz ^ 0), 



^E 



exp(/3£n2) - z' 



(5) 



f^ra is the number of particles in the kinetic ground state 
but in excited magnetic states. 



^E 



exp {fBion) — z ' 



(6) 



and Nmk is the number of particles in excited states of 
both magnetic and kinetic degrees of freedom. 



II. STRUCTURE OF THE LEVEL 
POPULATION 

Let us consider a charged gas of bosonic particles con- 
fined in a 3D box of volume AL, where A is the area 
in the x~y plane and L is the length in the z direction. 
Away from the edges of the box, the energy of a charged 
particle in a uniform ma gne tic field B pointing in the z 
direction is quantized as [30[ 



En{nz) = hw\n+ - ] + 



en. 



(1) 



where lo = qB/mc is the cyclotron frequency, e = 
H'^tt'^ /2mL'^ characterizes the level spacing in the z di- 



N^k = dJ2 E 



■'^ exp {/3enl + ^oju) 
The term A^^ can be approximated by 
A^fc ~ dr]kgi/2{z)- 



(7) 



(8) 



where 7?,. = (7r/4)i/2(/3e)-i/2 and g^{z) = ^Vj" 
stands for the polylogarithm function of order n [3l| . 
Here, the identity hmA^o J2jLi ^fi^j) = /o°° fi^) has 
been used assuming that A = \/pe is sufficiently small. 
Similarly, Nm and N^k can be approximated as: 



Nm ^ drirngiiz) 



(9) 
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and 

Nmk - dT]mVk93/2{z), (10) 

where ijm = (Z?'^)^^- We can group these approximations 
to write the total particle number as 

N 

-J = 9o{z) + Vkgi/2{z) + Vmgi{z) + VkVm93/2{z)- (n) 

As noted by the authors of Ref. in a similar eon- 
text (with harmonic instead of hard- wall confinement), 
the expression Eq. ([TT|) is an excellent and controlled 
approximation for the whole temperature regime. We 
note that, in this equation, the second, third, and fourth 
terms become dominant in the respective cases rjk 3> 1, 
rim > 1 and both 77^, ri,n S> 1. 

The factor d appearing in all the occupation numbers, 
Eqs. (|4|-(fTT|). reveals the existence of Landau degen- 
eracy, which in particular reflects a multiplicity of one- 
particle ground states where bosons are in both the ki- 
netic and magnetic ground state {n = = 0). 

The second term, Nk, represents the occupation of par- 
ticles lying in the ground Landau manifold but in ex- 
cited kinetic states (with motion in the z direction). We 
point out that this term has often been neglected in stud- 
ies of the CBG. We show below that Nk plays a most 
important role in the thermodynamic properties of the 
CBG under large magnetic fields. The third term, Nm, 
which accounts for those bosons which are in the kinetic 
ground state but in excited magnetic states, will not be 
studied here because it acquires importance only for very 
small magnetic fields, that is, when the energy of the first 
magnetic excited state is much lower than the energy of 
the first kinetic state (w <C s). This would describe a 
two-dimensional boson gas which will not be considered 
here. The last term, Nmk, corresponds to the conven- 
tional term in which the bosons are in both magnetic 
and kinetic excited states. It is the dominant term at 
high enough temperature, where the 3D limit is recov- 
ered {f]m,Vk > 1). 

Finally, we note that the order of the polylogarithm 
functions appearing in Eq. (|lip . 5^/2(^)1 corresponds 
to the population of excited states of free particles in d 
dimensions. 



III. NUMERICAL RESULTS FOR iVo AND A^fe 

For zero magnetic field and vanishing e the CBG re- 
sembles a 3D free Bose gas, showing BEC at the usual 
critical temperature. In this case, a discontinuous phase 
transition occurs [3^ . However, for large values of the 
ma gnet ic field (w ^ e), the transition becomes diffuse 
ITU and only crossover temperatures can be defined. 
Specifically, due to the presence of the Nk term, two dif- 
ferent crossover temperatures will be identified. 

It is quite remarkable that, for these high magnetic 
fields, Nk gives a significant contribution to A'', as re- 
vealed in Fig. [1] where we show the relative weight of 




T 



Figure 1. (Color online) Computed relative weight of No for 
iV = lO'^ and lj = 10^ (black lines), to = lO'^ (red -light gray- 
lines) and io — 10* (blue -gray- lines), in units of e. Temper- 
ature is given in units of e/kB- The conventional approach 
(which ignores Nk) is plotted with dashed lines. The effect of 
including the Nk in the relative weight of No is plotted with 
solid hues. The limiting function {1 + rjk)~^ is represented 
with black dotted lines (it overlaps with the solid blue (gray) 
line) [see Eq. (fT^ J. 



No{T) described by both the conventional and extended 
approaches {Nk — and Nk ^ 0, respectively) for sev- 
eral values of the magnetic field. From this figure it can 
be noted not only that the phase transition becomes dif- 
fuse, but also that the macroscopic occupation of bosons 
with kinetic excited states lying in the ground Landau 
manifold is very important at large magnetic fields. 

We note that, within this limit of large magnetic fields 
and for z ^ 1 (i.e. /3|/x| = —fipi ^ 1, high tempera- 
ture limit for the effective ID system), the temperature 
dependence of Nq/N tends to 



as shown in Fig. [TJ In fact, this limit function cannot be 
distinguished from the case a; = 10^ depicted in Fig. [TJ 

The origin of the discrepancy shown between the con- 
ventional and the extended approach (which increases 
with the magnetic field) can be traced to the relative 
weights of the different terms entering Eq. (fTTj). which 
explicitly includes the contribution of Nk- Take for ex- 
ample uj = 10'' and N = 10^. In this case, as shown in 
Fig. [21 Nk represents a significant (and occasionally dom- 
inant) amount of the total particle number. Moreover, 
two different crossover temperatures appear when either 
Nk or A^o begin to grow substantially. In particular, the 
phase transition that occurs at a lower temperature re- 
mains somewhat more diffuse than the high-temperature 
one, because of its effective ID nature. 
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Figure 2. Computed relative weights of A'o (solid line), A'^^ 
(dashed line) and N^k (dotted line) for lo = 10^ and A'' = 10^ . 
See text for details. 



IV. SPECIFIC HEAT 



Figure 3. (Color online) Specific heat at constant volume as a 
function of temperature. The extended and conventional ap- 
proaches are plotted with solid and dashed lines, respectively. 
We have taken = lO'' and uj = 10^ (black lines), uj = 10^ 
(red -light gray- lines) and a; = 10* (blue -gray- lines). See 
text for details. 



The essential features of BEC as a phase transition are 
clearly exhibited in the behavior of the specific heat. In 
particular, we consider the specific heat at constant vol- 
ume, defined by Cy = — Z?'^ {dU/dl3)y. In the extended 
approach, the internal energy can be expressed as 

U{(3,z) = -,^fe/3-ig3/2(z)-f- —77^,7^/3-155/2(2). (13) 



When no magnetic field is present, Cy tends to the 3D 
free Bose gas, displaying a sharp maximum and a dis- 
continuous first derivative. The effect of a nonzero mag - 
netic field is that of smoothing out that maximum [14 1. 
This can be clearly appreciated in Fig. |3l Although the 
phase transition remains diffuse in the presence of a mag- 
netic field, as previously discussed, one can still associate 
a critical temperature to the maximum of the specific 
heat, as in the free gas case IJ]. In Fig. |3] it can also 
be seen that, for uj = 10^, no significant deviations from 
the conventional approach are found except for one: the 
inclusion of the Nk terms tends to round off the rela- 
tively sharp maximum shown in the specific heat calcu- 
lated from the conventional approach. On the contrary, 
for u! = 10^, significant differences between the conven- 
tional (dashed line) and extended (solid line) persist in a 
wide range of temperatures. Specifically, for large mag- 
netic fields (w ^ e) and high temperatures (/3a; <C 1, 
or ry,„ ^ 1), Cy/N approaches the classical value of 3/2. 
Within this regime, the classical limit can be obtained by 



taking N ~ drjkrji^z and U ~ ^drjkrjaiP~ z, which results 
in [/ = 3N/2(3. This is the classical limit for a system in 
three dimensions. 

It is interesting to note that signatures of a double 
phase transition appear in the specific heat. In particu- 
lar. Fig. [3] shows a plateau-like structure in the specific 



heat when the effective number of particles, N/d, is of the 
order of 10"^ (solid blue -gray- line). This plateau spans 
the temperature range e <^ ksT <^ uj. As T decreases, 
the first phase transition, observed at a high tempera- 
ture, reveals the transition from states which have both 
kinetic and magnetic excitations to those which are in 
excited kinetic states but in the ground Landau mani- 
fold. The second phase transition shows how the abso- 
lute ground state (kinetic and magnetic) is reached. In 
fact, the existence of two crossovers can be understood 
from Fig. [21 This is in accordance with the large slope 
of the Cv{T) curve observed for the high-temperature 
crossover, compared to the low-temperature one. On the 
contrary, the specific heat calculated within the conven- 
tional approach (that which ignores Nk) does not show 
signatures of double phase transitions, as expected. 

We remark that, in the plateau e ksT <^ ut, a. 
constant value for the specific heat is obtained between 
the two crossover temperatures. Within this temperature 
range, N ~ drj^z and U ~ ^drj^P^^ z. Finally, we get 
U ~ N/2(3. Thus, the plateau of the specific heat lies 
in the value Cy = Nks/'i, as shown in Fig. [31 which 
corresponds to that of a Bose gas in one dimension. We 
note that, although the value of the specific heat in the 
intermediate temperature plateau is independent of the 
degeneracy d, the details of the low and high temperature 
behavior depend on the particular value of d. In fact, 
we recall that the very existence of the plateau relies 
essentially on the degeneracy being much greater than 
unity. 
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V. CONCLUSIONS 

In this work we have revisited the free charged Bose gas 
focusing on the case of a large magnetic field. Although 
a sharp Bose-Einstein condensation does not occur, two 
diffuse phase transitions (or crossovers) are found. As the 
temperature is decreased, the first crossover reflects the 
enhanced occupation of states in the lowest Landau level. 
Then a range of intermediate temperature is found where 
the system behaves has an effectively one-dimensional 



system in its high-temperature regime. At even lower 
temperatures, the system undergoes a ID Bose conden- 
sation, which is known to be a diffuse phase transition. 
To account for this two-step condensation, it has proved 
essential to include the occupation of states in the low- 
est Landau level but with arbitrary kinetic energy in the 
direction parallel to the magnetic field. 

This work has been supported by MINECO (Spain) 
through grants FIS2010-21372, CTQ2008-02578/BQU 
and the Juan de la Cierva program (PB). 
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